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Abstract 

Naturally occuring or man-made systems displaying periodic spatial modulations 
of their properties on a nanoscale constitute superlattices. Such modulated struc- 
tures are important both as prototypes of simple nanotechnological devices and as 
particular examples of emerging spatial inhomogeneity in interacting many-electron 
systems. Here we investigate the effect different types of modulation of the system 
parameters have on the ground-state energy and the charge-density distribution of 
the system. The superlattices are described by the inhomogeneous attractive Hub- 
bard model, and the calculations are performed by density-functional and density- 
matrix renormalization group techniques. We find that modulations in local electric 
potentials are much more effective in shaping the system's properties than modula- 
tions in the attractive on-site interaction. This is the same conclusions we previously 
(Phys. Rev. B 71, 125130) obtained for repulsive interactions, suggesting that it is 
not an artifact of a specific state, but a general property of modulated structures. 
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1 Introduction 



Superlattices are characterized by a periodic modulation of system properties 
on a length scale that is larger than the lattice constant of the material they 
are made of. Superlattices are of great technological interest because the mod- 
ulation can be used to design the system's optical, electrical, magnetic and 
transport properties [1 .2,3,4, 5j . Since the modulation typically takes place on 
a length scale of nanometers, superlattices also provide a particular realization 
of systems with nanoscale spatial inhomogeneity, which occurs spontaneously 
in many strongly correlated systems [5117115] . 

A simple, but nontrivial, model of a superlattice is the inhomogeneous one- 
dimensional Hubbard model (1DHM) for spin 1/2 fermions on a lattice, with 
Hamiltonian 

i,tr i itr 

where either the on-site interaction C/j, or the on-site potential V{ or both are 
periodically modulated. Both for designing artificial superlattices and for the 
study of spontaneously emerging periodic structures, it is important to know 
which of these two types of modulations is more effective in shaping the system 
properties. 

In Ref. [H] we investigated this question for repulsive interactions, and found 
that the properties of the superlattice are much more sensitive to local changes 
in the electric potential than to local changes in the strength of the interpar- 
ticle interaction. This is a surprising result, as intuitively one would expect 
that in a strongly correlated system, such as described by the 1DHM, inter- 
action effects should dominate. Indeed, in most applications of the Hubbard 
model (not only to superlattices, but also, e.g., to HTSC cuprates) the on-site 
potential is not even considered, or simply assumed to be spatially constant 
and absorbed in the chemical potential. In the present paper we extend this 
investigation to attractive interactions, i.e., U < 0. The negative U Hubbard 
model has been used to describe superconductivity of strongly correlated elec- 
trons, but the same model also describes attractive interactions between cold 
atoms in optical traps [10J. 



2 Comparison of modulated structures 

Many different modulation patterns exist, and many different ways to compare 
them are conceivable. What we require for our present purposes are quanti- 
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tative criteria measuring the relative impact of modulations in the on-site 
potential, as compared to modulations in the on-site interaction. 

In Ref. [9] we introduced two such criteria, which are both based on separate 
calculations for a 1DHM with no modulations (a spatially uniform system) 
and for a 1DHM in which Ui and Vi are modulated with same periodicity but 
different amplitudes. In the present paper, dedicated to attractive interactions, 
we choose a periodically repeated pattern of Lu sites with interaction [/, < 
and potential Vi, followed by L Q sites with Ui = Vi = 0. The modulations of 
Ui and of Vi have distinct effects on the system properties. The ratio of the 
modulation amplitudes at which the two effects cancel each other, returning 
effectively to the properties of the uniform system, provides a quantitative 
measure of their relative impact on the system properties. 

Two important ground-state properties of a modulated structure are its total 
energy and its density profile. These quantities form the basis for the two 
criteria we employ. In the energy criterium, we vary the on-site potential until 
the doubly modulated lattice has the same ground-state energy found in the 
uniform system. The density criterium consists in finding that potential that 
produces the same density profile obtained in the open uniform system. Exact 
cancellations are harder to obtain for the density than for the energy, but 
cancellations to within the error bars of our numerical methods are all that is 
required for assessing the relative impact of modulating i> f and Ui. 

We stress that we do not expect such cancellations to occur frequently in 
nature. We simply use them here as a convenient way to quantify the relative 
impact of both types of modulation. Essentially, we are simulating numerically 
the sensitivity of the energy and the density to changes in the interaction and 
the potential, i.e., sampling certain response functions. 



3 Computational methods 

We calculate the energies and densities by means of two independently de- 
veloped and implemented numerical techniques: the density-matrix renormal- 
ization group (DMRG) [TT][T2] and density-functional theory (DFT) within 
the Bethe-Ansatz local-density approximation (BA-LDA) p3l!15lll6p!7] . Both 
methods provide energies and density profiles of the Hamiltonian (Tj[|), but 
are complementary in scope. DMRG provides high-precision results, but at 
considerable computational expense, in particular for systems with reduced 
spatial symmetries or periodic boundary conditions. BA-LDA easily handles 
systems of almost arbitrary size and its performance does not significantly 
depend on presence or absence of spatial symmetries or choice of boundary 
conditions. Its intrinsic accuracy is, however, limited to a few percent (with 
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Table 1 

Ground-state energy of superlattices with open boundary conditions, L sites and 
./V electrons, subjected to spatially modulated on-site interaction {J, and spatially 
constant on-site potential Vi = 0, obtained with DMRG and with DFT/BA-LDA. 
The modulation pattern alternates Ljj sites of interaction U < with Lq sites with 
interaction {7 = 0, simulating a superlattice composed of alternating superconduct- 
ing and normal metals. The last column is the deviation of the BA-LDA energies 
from the DMRG ones (in percent). 
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Lu 
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E BMRG /t 


E* A - LDA /t 
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30 


-3 


3 


2 


-69.473 


-70.379 
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-337.92 


-340.88 


0.88 


300 


150 


-6 


10 


10 


-515.18 


-532.79 


3.42 


300 


300 


-3 


10 


10 


-556.07 


-560.41 


0.78 



energies typically being slightly more accurate than densities) and it is not 
easy to consistently improve on these values. 



Within density-functional theory, the ground-state energy and density dis- 
tribution are obtained by self-consistent solution of a set of single-particle 
equations, the Kohn-Sham equations, in the presence of effective potentials 
constructed such as to recover the density distribution of the interacting sys- 
tem. The application of DFT to the Hubbard model was initiated in Refs. 
[T3] . Recent progress is due to the development of a simple and accurate 
Bethe-Ansatz based parametrization for the model's correlation energy in the 
uniform limit, which can be used to construct a local-density approximation 
(LDA) [T4"lll5l)16j . This parametrization was further developed into a fully 
numerical scheme [17] . which has been applied both to repulsive [T7j and at- 
tractive [10J interactions, in the context of ultracold fermionic atoms confined 
in optical traps. In the same context, a simple analytical parametrization of 
the LDA for attractive interaction was also proposed [18] , but applied only in 
conjunction with a similar LDA for the kinetic energy, leading to a Thomas- 
Fermi-like scheme |18|. 



Here we employ the parametrization for the correlation energy used in Ref. 
|18j . in conjunction with the Kohn-Sham approach of Refs. [T4~|15|ll6j . to ob- 
tain many-body energies and densities. In these calculations a local-density 
approximation is made only for the correlation energy. Previous work [10] . as 
well as our present results, indicate that the LDA (in either its fully numerical 
or parametrized forms) for attractive interactions continues to produce ener- 
gies within a few percent of more precise DMRG data. This is illustrated in 
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Fig. 1. Amplitude v of the modulation of the on-site potential that produces the 
same ground-state energy in a superlattice whose interaction is already modulated 
with amplitude U, as obtained in an unmodulated system with interaction U on all 
sites and v un if = 0. Full squares: DMRG data. Plus: DFT/BA-LDA data. System 
parameters: L = 60 sites, N = 30 fermions, Ljj = 4 interacting sites, alternating 
with Lq = 6 noninter acting sites. 

Table [TJ DMRG densities are less well reproduced by LDA densities at attrac- 
tive than at repulsive interactions, but Figs. Q] and [2] show that the differences 
are far too small to invalidate our conclusions regarding the effect of different 
modulation patterns. 

In view of the complementary aspects of DFT and BA-LDA , we have, both 
in Ref. [9] and in the present paper, used BA-LDA to efficiently explore large 
regions of parameter space. In separate calculations we have then used DMRG 
to fine-tune parameters in the vicinity of the values provided by BA-LDA, and 
to confirm results with more significant digits. The DMRG calculations were 
performed using the finite system algorithm under open boundary conditions 
[TT] with truncation errors kept of the order of 10~ 6 or smaller. 



4 Attractive Hubbard model: results 



Figure [T] displays representative results from the energy criterium. In these 
calculations we have chosen a fixed modulation of the interaction, alternating 
'layers' of Ljj = 4 sites with Ui = U and 'layers' of L = 6 sites with Ui = 0, 
and changed the amplitude of the modulation of Vi until the doubly modulated 
system has the same energy as the unmodulated system (with interaction U 
on all sites and v uni f = 0). If a different constant value of v uni f is chosen in 
the uniform system, the compensating values of Vi in the doubly modulated 
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Fig. 2. Amplitude v of the modulation of the on-site potential that produces the 
same ground-state density distribution in a superlattice whose interaction is already 
modulated with amplitude U, as obtained in an unmodulated system with interac- 
tion U on all sites and v un if = 0. Full squares: DMRG data. Plus: DFT/BA-LDA 
data. System parameters: L = 160 sites, TV = 80 fermions, Ljj = 6 interacting sites, 
alternating with Lq = 10 noninteracting sites. 

system are simply shifted up by that value of v uni f, without changing our 
conclusions. 

Figure [1] shows that for physically realistic values of U the effect of a mod- 
ulation of U on the ground-state energy is reproduced by a modulation of v 
that is between 10~ 2 and 2 x 10 _1 times smaller. Clearly, the ground-state 
energy depends much more sensitively on modulations of the on-site potential 
than on modulations of the on-site interaction. Qualitatively, this is the same 
conclusion obtained previously for repulsive interactions 0. Quantitatively, 
the effect of attractive interactions is even weaker than that of repulsive inter- 
actions, as can be seen by comparing the present Fig. [1] with the upper panel 
of Fig. 4 of Ref. [9]: an even weaker modulation of the potential is sufficient 
to change the energy by the same margin as a modulation of the interaction. 

In Fig. [2]we turn to the density criterium, i.e., investigate for what amplitude 
of modulation of an repulsive on-site potential Vi it cancels the effect of a 
modulation of the attractive on-site interaction [/, on the density profile. The 
density is more sensitive to small variations in the system parameters, and 
in general it is more difficult to satisfy the density criterium than the energy 
criterium. In practice, we have performed a large number of BA-LDA calcu- 
lations within DFT, until finding values of Vi that produce approximately the 
same density profile. In the vicinity of these BA-LDA values we have then 
performed a few DMRG calculations, for selected values of U, to confirm the 
observed trend. This procedure is motivated by the large difference in comput- 
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ing time required by both methods (seconds to minutes, as compared to hours 
or days, depending on the system parameters and the required precision). 

The data in Fig. [2] show the same trend as those in Fig. [U indicating that 
the density profile of the superlattice is about one order of magnitude more 
sensitive to changes in the potential than in the interaction, even for quite 
large values of U. This is the same finding obtained in Ref. [9] for repulsive 
interactions. Since the nature of the ground state at attractive and repulsive 
interactions is very different, we infer that we are observing a general feature 
of periodically modulated systems, not a specific consequence of one type of 
ground state. The ratio between compensating potentials and interactions is, 
however, not the same for the energy criterium and for the density criterium, 
and also differs, for each criterium, for attractive and for repulsive interactions, 
which indicates that there is no universal energy scale associated with the 
compensation. 



5 Conclusions 



Density-functional and density-matrix renormalization group calculations have 
been performed for superlattices with different modulation patterns, in the 
presence of attractive interactions. Results for ground-state energies and charge- 
density profiles, obtained for different lattices and with DFT or DMRG, are 
mutually consistent and all point towards the same conclusion: modulations 
in the on-site potential are much more efficient in shaping the system prop- 
erties than modulations in the on-site interaction. This observation has two 
distinct consequences, depending on whether ones interest lies in the design 
of materials and devices, or in the understanding of spontaneously emerging 
spatial inhomogeneity. 

If one wants to design a superlattice with specific electronic properties, e.g., 
as part of a nanoscale device, the key parameter to modulate is the local elec- 
tric potential. Spatial variations of the effective particle-particle interactions 
are of secondary importance. If, on the other hand, one is given a system 
in which nanoscale spatial inhomogeneity appears spontaneously, such as in 
many HTSC cuprates, our results imply that the analysis of origins and con- 
sequences of the modulation in terms of strong-correlation (large U) effects 
alone is incomplete, and must necesssarily be accompanied by accounting for 
local electric fields. 
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